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Phase diagram of microcanonical ensembles of self-attracting particles is studied for two types of 
short-range potential regularizations: self-gravitating fermions and classical particles interacting via 
attractive soft — (r 2 + r 2 ,) -1 ^ 2 Coulomb potential. When the range of regularization is sufficiently 
short, the self-attracting systems exhibit gravitational or collapse-like transition. As the fermionic 
degeneracy or the softness radius increases, the gravitational phase transition crosses over to a 
normal first-order phase transition, becomes second-order at a critical point, and finally disappears. 
Applicability of a commonly used saddle-point or mean-field approximation and importance of 
metastable states is discussed. 
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I. INTRODUCTION 

Ensembles of particles interacting via a long-range non- 
integrable attractive potential, U(r) = Ar~ a , A < 0, 
and < a < 3, are known to exhibit gravitational 
phase transition between a relatively uniform high energy 
state and a low-energy state with a core-halo structure 
§1, II- It has also been 

established |4|, ]|, || ]7], g, |l(J that if the interaction 
potential has some form of short range cutoff, the den- 
sity and all other physical quantities that characterize 
the core-halo state are finite, while for a bare r~ a at- 
tractive potential the collapse results in a singular state 
with an infinite density, entropy and free energy. The 
gravitational phase transitions are known to exist in mi- 
crocanonical, canonical, and grand canonical ensembles; 
yet the details of the phase transition and the structure 
of the core-halo state are ensemble- dependent. In this 
paper we shall essentially consider microcanonical ensem- 
bles; for long-range interacting systems they are known 
to have richer phenomenology and allow for states (such 
as ones with a negative specific heat) that are inaccessi- 
ble in both canonical and grand canonical ensembles. In 
addition, the microcanonical ensemble is the most funda- 
mental since the notion of thermostat is ambiguous for 
long-range interacting systems. 

A typical entropy vs. energy plot, describing a grav- 
itational phase transition in a microcanonical ensemble, 
is shown in Fig. [l], lower plot ||, ||, [?], ||, [l0|] . High-energy 
and low-energy branches correspond to the uniform and 
to the core-halo states, respectively; their intersection E* 
marks the point of the phase transition. Metastable uni- 
form and core-halo states, shown by dashed lines, exist 
in the energy intervals [E\,E*] and [E* , E2], correspond- 
ingly. When the energy is larger than E2, the system 
can be only in the stable uniform state. The system still 
stays in this uniform state after energy is decreased past 
E* and the uniform state becomes metastable. When the 
energy reaches E\, the uniform state becomes unstable 
and the system undergoes a collapse to a stable core-halo 
state. During such a collapse the entropy is discontin- 
uously increased, and the macroscopic rearrangement of 
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FIG. 1: Sketch of an entropy vs. energy plot for a mi- 
crocanonical self-attracting system with a short-range cutoff 
which exhibits a gravitational (G) phase transition. Stable 
states and metastable states are shown by solid and dashed 
lines, respectively. Points E\ , E* , and E2 denote the collapse 
energy for the uniform state, the energy of the phase transi- 
tion, and the energy for which the metastable core-halo state 
becomes unstable, correspondingly. A typical entropy vs. en- 
ergy plot for a microcanonical long-range interacting system 
with a normal (iV) first-order phase transition is also shown. 



the density profile occurs. Similarly, the core-halo state is 
stable below E* , metastable between E* and E2 , and un- 
dergoes a discontinuous transition reverse to the collapse, 
at i?2- This transition is sometimes called an explosion 
||, since it transforms the dense core into a relatively 
uniform mass distribution. 

In some sense, the way the gravitational phase transi- 
tion occurs in a microcanonical ensemble resembles a hys- 
teresis phenomenon that takes place during a first-order 
phase transition in a canonical ensemble. For this reason, 
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it is sometimes called a gravitational first order phase 
transition. Yet in a microcanonical ensemble a normal 
first-order phase transition occurs without the hysteresis 
and metastable states. A distinct feature of a normal 
microcanonical first-order phase transition in a small or 
long-range interacting system is a convex dip in other- 
wise concave continuous entropy plot (see, for example, 
U, H|, [[§) sketched in Fig. [y, upper plot (N). Differ- 
ences between gravitational and normal first-order phase 
transitions are evident immediately: for a gravitational 
phase transition the microcanonical inverse temperature 
/3 = dS/dE is discontinuous at E* while for the normal 
first-order phase transition dS/dE is always continuous. 
In addition, in a gravitational phase transition, the uni- 
form and the core-halo phases cannot coexist, while for 
the normal first-order phase transition the phases do co- 
exist in a range of energies corresponding to the convex 
dip. 

It has been noticed || ||, [?J ||, [l(J that the gravita- 
tional phase transitions in self-attracting systems exist 
only when the range of the cutoff is sufficiently small. 
When the effective cutoff radius is increased, the range of 
existence of metastable states [E\, E2] shrinks and finally 
disappears along with the dS/dE discontinuity at E*. 
The resulting entropy vs. energy plot becomes continu- 
ous and qualitatively resembles a corresponding plot for a 
system with a normal first-order phase transition |^|, [l0| . 
This was observed in self- gravitating systems with vari- 
ous short-range regularizations: central excluded volume 
JjJ , hard sphere repulsion for individual particles || , soft- 
core interaction potential || , and exclusion due to Fermi- 
Dirac statistics [|, [To| . 

In this paper we study how the behavior of self- 
attracting systems depends on the effective short-range 
cutoff radius. We find that once the cutoff is increased 
above a certain value, the gravitational phase transition 
crosses over to a normal first-order phase transition, char- 
acterized by a convex dip in the entropy vs. energy plot 
and associated with an energy interval with a negative 
specific heat. As the effective cutoff radius is increased 
even further, the system reaches a critical point where the 
first-order phase transition is replaced by a second-order 
one, and for even larger cutoff radii there is no phase tran- 
sition at all. To reveal that such a phase diagram is a 
generic feature of all self-attracting systems, we consider 
two different examples: an ensemble of self-gravitating 
particles with Fermi-Dirac exclusion statistics, and a par- 
ticle system with classical statistics interacting via a soft 
Coulomb potential — (r 2 + Tq) -1 / 2 . Since both of these 
examples have been studied in the past, in the following 
two sections we present only short outlines of the deriva- 
tion of the main formulas and refer readers to || ||, [| [T(J 
for a more detailed description. After analyzing these 
two examples, we discuss the validity of the saddle-point 
or mean-field approximation, which is commonly used to 
study self-gravitating systems; a conclusion which sum- 
marizes the presented results completes the paper. 



II. SELF-GRAVITATING FERMIONS 

We consider a microcanonical ensemble of N ^> I iden- 
tical unit mass fermions confined to a spherical container 
of radius R. A simple combinatorial analysis gives for the 
number of microscopic states W of a Fermi-Dirac system 
with a phase space distribution function /(r, p): 

W\f\K 

[nlnn+ (1 - n) ln(l - n)]d 3 rd 3 p > .(I) 



exp 



(2nh) 3 



Here n(r, p) = (27rfi,) 3 /(r, p)/g; g is a number of internal 
degrees of freedom, g — 2s + I for a particle with spin 
s. The logarithm of Eq. is the Fermi-Dirac entropy. 
The most probable distribution expected at equilibrium 
is obtained by maximizing Eq. (Q) with respect to n with 
the constrained total energy E and number of particles 
N. The validity of this mean-field approach is discussed 
in Sec. 
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Expressing the particle density through the 
Laplacian of the gravitational field by using the Poisson 
equation, it is found that the mean-field gravitational 
potential </>(r) satisfies an equation of the form 



A0(r) 



4nGg f d 3 p 

{2-Khf J l + e lP 2 /2+4>(r)-v]/T> 



(2) 



where G is the gravitational constant. The temperature 
T(E, M) and the chemical potential v(E, M) appear as 
Lagrange multipliers associated with the conservation of 
energy and mass. We assume a spherical symmetry of 
the mean-field potential, introduce a dimensionless ra- 
dius £ = r{Gg^/WT fitti 3 ) 1 / 2 , a non-negative dimension- 
less potential = [4 > { r ) ~ < fi(fy]/T, and a positive con- 
stant k = exp{[0(O) - v]/T}. Then Eq. (|) is reduced to 
the form Q: 
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Here Ii/ 2 denotes the Fermi integral 
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Boundary conditions if>(0) = and ^'(0) = follow from 
the definition of the potential if)(t;) and spherical sym- 
metry, respectively. Using the Gauss theorem for the 
gravitational potential at the boundary of the spherical 
container R, the number of particles constraint can be 
expressed as 



dr 



GM 
TR2' 



(5) 



Introducing a dimensionless maximum radius a = 
R{GgV8T/Trfi 3 ) 1/2 , we reduce Eq. (|) to the form 

«V(«) = M 2 , (6) 
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where fi = y / 8g 2 G 3 MR 3 /ir 2 h 6 is a degeneracy param- 
eter H , which is proportional to the ratio of the grav- 
itational Fermi energy (defined in Eq. (]?]) below) to 
the average gravitational energy GM 2 /R to the power 
3/2. In the large-/i limit the system becomes completely 
non-degenerate or classical, and Eq. (||) then reduces 
to the Boltzmann-Poisson equation. Yet for any finite \i 
and sufficiently low energy the system is degenerate and 
mimics the structure of a white dwarf star In the 

completely degenerate limit T — ► 0, the equation for the 
gravitational potential becomes the Lanc-Emden equa- 
tion for a polytrope of index n = 3/2 The density 
goes to zero at a finite radius i?» related to the total mass 
M by the relation MRl = 91.869fi 6 G*- 3 (2/ 5 ) 2 . The en- 
ergy of this fermion ball is given by the Ritter's formula 
E = — (3/7)GAf 2 /i?». Hence the minimal possible en- 
ergy of a self-gravitating degenerate gas, i.e. the gravi- 
tational Fermi energy is 



= -5.98 icrV /3 



G 2 M 7 / 3 

h 2 ■ 



(7) 



Note that contrary to the completely degenerate (E — 
E m i„) case, the density of the self-gravitating Fermi gas 
for E > E min decreases like r" 2 at large distances, 
similarly to the classical self-gravitating isothermal gas. 
Hence the container is necessary to confine the system 

As in |6|], we introduce a dimensionless energy e, a 
dimensionless inverse temperature f3, and the entropy per 
particle s, and express them through a, [a, tp and k as 
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GM 2 
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-e/3 + tp(a) + /3 + lnfc 



2a 6 
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Hence, all the thermodynamics parameters e, s, and (3 
are single- valued functions of the degeneracy parameter 
fi and the uniformization variable k, which continuously 
parametrizes the functions V'(C)- We select a degeneracy 
(i, and for each k we numerically integrate Eq. (||) for 
< £, < a, where a is determined by Eq. (g). Once 
?/>(£) is computed, e(k), f3(k), and s(k) are calculated 
using Eq. (|S|). The curves giving the entropy s(k) vs. 
energy e(k) and the inverse temperature (3{k) — ds/de 
vs. energy e(k) are thus defined in a parametric form for 
various values of the degeneracy parameter fj, (see Figs. 

In the classical fj, — > +oo limit, the (e,/?) plot forms a 
spiral which winds up indefinitely around a limit point 
corresponding to the singular isothermal sphere with the 
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FIG. 2: Plots of entropy per particle s(e) (upper plot) and 
dimensionless inverse temperature /3(e) = ds/de (lower plot) 
for fi = 10 9 . A change of stability of the corresponding state 
occurs each time the (e, f3) spiral has a vertical tangent and 
the (e, s) plot has a cusp. A mode of stability is lost at the 
vertical tangent if the curve rotates anticlockwise and gained 
if the curve rotates clockwise jL8|. Therefore, the branches 
(1) and (7) are stable, (2) and (6) are unstable against one 
mode, (3) and (5) against two modes and (4) against three 
modes (branch (7) is not represented in the figure but is the 
continuation of branch (6) after the turning point of energy at 
e ~ 147) . There exists values of energy at which the branches 
with the same degree of instability have the same entropy. 
It occurs at the crossing points of the s vs e plot and cor- 
responding vertical lines in the /3 vs e plot. As /i — * +oo, 
there are more and more crossing points at energies e n con- 
verging to the value e = 1/4 of the singular isothermal sphere. 
These points (e„,s n ) can be associated with points of phase 
transitions; however, these "phase transitions" occur between 
unstable states and are therefore unphysical. 



energy e = — 1/4 and inverse temperature (3 = 2 |17| • For 
very high but finite values of fi, such as /i = 10 9 in Fig. 
the system is almost non-degenerate and its (3 vs. e plot 
looks very similar to this infinite winding spiral. Yet as 
fj, is finite, the spiral is finite as well and does unwind 
after a turning point. Each point with a vertical tangent 
(where dft/de = oo) on the [3 vs. e spiral corresponds to a 
cusp in the s vs. e plot. However, all the exotic features 
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FIG. 3: Plots of entropy per particle s(e) — 5 (top) and 
dimensionless inverse temperature /3(e) = ds/de (bottom) for 
fj, — 10 4 . Unstable state is shown by dashed lines. 



FIG. 4: Plots of inverse temperature /3(e) = ds/de for (left to 
right) /x = 2.65 x10 s , 2.67 xlO 3 , 2.7 x 10 3 , and 2.75 x 10 3 . The 
gravitational phase transition disappears for /i < [i gr — 2.67 x 
10 3 and is replaced by a normal first order phase transition. 



of the high-/x plots, such as multiply-winding spirals and 
corresponding multiple kinks are related to unstable and 
therefore physically unrealizable states |jl8| . 

As the degeneracy parameter /z is decreased, the num- 
ber of unstable states is decreased as well. As illustrated 
in Fig. H[ for \i = 10 4 there is only one unstable state, 
shown on the plot by a dashed line. Yet for /i = 10 4 
the system is still highly non-degenerate and exhibits 
all signs of the gravitational phase transition. The low- 
energy branch of the plot corresponds to the core-halo 
state which exists for < k < k% and e m i n < e < £2, 
where e m in = —6.42 10 _2 ^( 2 / 3 , is the dimensionless grav- 
itational Fermi energy. Numerical evidence suggests that 
at the 62-cusp, both energy e(fc) and entropy s(k) behave 
as 0(k — fe) 2 , which explains the divergence of df3/de, 
evident in the plot. For k2 < k < k\ the entropy s(k) 
is not even a local maximum || |hJ so the correspond- 
ing state, marked by the dashed line in Fig. ||, is un- 
stable. When the curve approaches the second cusp at 
(-\ = e(fci), both energy e and entropy s go through a 
0[(k — ki) 2 ] asymptotics again. For k > k\ the equilib- 
rium states, now belonging to the uniform phase, are at 
least locally stable. 

When [i gets smaller and the degeneracy becomes more 
important, the interval between t\ and €2 decreases and 
finally disappears. In Fig. || we plot /3 vs. e and seek 
/x such that the (3(e) curve loses its characteristic for 
large /i backbend. This can be viewed as a complete 
straightening of the classical spiral and happens when 
/i = jj, gr ps 2.67 x 10 3 . For /1 < fi gr the system ex- 
hibits a normal first-order phase transition. To illustrate 
this, in Figs. ^-|^ we present entropy and inverse tempera- 
ture plots for /1 = 10 3 ; signs of microcanonical hrst-order 
phase transition such as convex entropy dip and region 
of negative specific heat where d 2 s/de 2 > are clearly 



visible. 

As we decrease the degeneracy parameter fj, even fur- 
ther, the convex dip in the entropy vs. energy plot and 
correspondingly the interval where d 2 s/de 2 > get nar- 
rower, and for [i = fi cr « 83 they disappear (Fig. 0). At 
M = Mcr the equation d 2 s/de 2 = has only one real root, 
e = e cr w —0.5, a critical point where the two phases 
become indistinguishable and the heat capacity diverges. 
This corresponds to the line of first-order phase transi- 
tions in (e, fj,) space, terminated by the critical point at 
(e cr ,Ucr) where the phase transition is second order (see 
Fig. 0)- 

For [i < fi cr , the inverse temperature /3(e) is a mono- 
tonically decreasing function of energy e and the sys- 
tem does not exhibit any phase transition. Therefore, 
as the degeneracy parameter decreases, the microcanon- 
ical ensemble of self-gravitating fermions consecutively 
exhibits gravitational, first-order, second-order, and no 
phase transition at all. Additional discussion of phase 
transitions in the self-gravitating Fermi gas can be found 
inRcf. 0. 



III. SELF-ATTRACTING PARTICLES WITH 
SOFT COULOMB POTENTIAL 

In this section, we consider phase transitions in another 
self- attracting system where the short-range cutoff is ex- 
plicitly present in the interaction potential. As in Sec. O, 
we consider a microcanonical ensemble of N 1 identi- 
cal unit mass particles confined to a spherical container of 
radius R, but now the particles obey classical statistics 
and interact via the attracting soft Coulomb potential 
—G/(r 2 + r 2 .) -1 / 2 . This potential has been used in vari- 
ous numerical simulations of self-gravitating systems and 
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FIG. 5: Plot of entropy s vs. energy e for jj, = 10 . The 
entropy presents a convex intruder between e and e . For 
an extensive system, this convex intruder is forbidden be- 
cause the system with energy e, e < e < e'" , would gain 
entropy by splitting in two phases with energy e and e'". In- 
deed, s(e = ae' + (1 - a)e"') < as(e') + (1 - a)s(e"'), where 
< a < 1 parametrizes the energy e of the systems in the 
phase coexistence range [e , e'"]. However, for a non-extensive 
system, such as a gravitational system, this argument does 
not hold and a convex intruder for the entropy is allowed in 
the microcanonical ensemble |p.3i. 



is also called the Plummcr potential. Since such poten- 
tial does not satisfy the Poisson equation, the equation 
for the mean-field density or potential cannot be reduced 
to a differential equation similar to Eq. (^). In this sec- 
tion, we use the integral equation approach suggested in 
Ref. ||. The entropy S(E) of the system is defined as 
the logarithm of the density of states with the energy E, 



S(E) 



xS 



M 




M 

n 

k=l 

G 



d 3 p k d 3 r k 
(27rft) 3 
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Following the steps described in Refs. @, H], we inte- 
grate Eq. (JsJ) on momenta, express the remaining con- 
figurational integral through a functional integral over 
possible density profiles p(r), apply the saddle-point ap- 
proximation, and introduce the dimensionless coordinate 
x = r/R, soft potential radius xo = r /R, and energy per 
particle e = ER/GM 2 . The saddle-point density profile 



FIG. 6: Plot of inverse temperature /3(e) = ds/de vs. en- 
ergy e for fj, — 10 . The existence of negative specific heats 
C = —f3 2 de/d[3 < and the convex intruder for the entropy 
are the signal of a normal first order phase transition and of 
the inequivalence of statistical ensembles. Indeed, for nonex- 
tensive systems, the region of negative specific heat is allowed 
in the microcanonical ensemble while it is forbidden in the 
canonical ensemble and replaced by a sharp phase transition 
(horizontal plateau). The temperature of the transition /3 t _1 is 
determined by the crossing point in the free energy J = s — f3e 
vs inverse temperature (3 plot (dashed line). Alternatively, it 
can be obtained by performing a Maxwell construction in the 
/3 vs e diagram, noting that f**(J3 — /3 t )de = J'" — J' = 
(the areas of the shaded regions are — AJ = J" — J' and 
A J = J'" — J"). It is also given by the slope of the straight 
line s h uii{e) in Fig. | (s h uii = ae + b with a = f3 t and 
b — J' — J'") [P-3[] . During the canonical phase transition, 
a latent heat Ae = e'" — e is released. 



p s (x) satisfies the following integral equation 
p s (x) = p F[p s (.),x], 



F[Ps(-),x] = cxp 



jO s (xi)/9 s (x 2 ) 



6 2 / / . ; v . - v-.. i- j- 



P0 



V(xi - x 2 ) 5 
F[p s (.),x]d 3 x 



d 3 xid 3 x 2 



(10) 



Neglecting an energy-independent constant, the saddle- 
point entropy per particle s(e) = S(e)/M is expressed 
as 



fl(e) = ™ln/3- / Ps (x)ln[/,(x)]d 3 x, (11) 



where (3, introduced in Eq. (^0|), is the inverse dimen- 
sionless temperature /3(e) = ds/de. In order to solve 
Eqs. (|l0|), we assume spherical symmetry of p, integrate 
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FIG. 7: Plots of inverse temperature /3(e) = ds/de vs. energy 
e for (left to right) p = 100, 90, 85, 80, and 50. The normal 
first order phase transition disappears for p, < p cr m 83. At 
the critical value p = /i cr the system passes by a second order 
phase transition. 



Eqs. (10) over angular variables and introduce a map 



p i+ i(x) = ap F[ Pi {.),x] + (1 - a) Pi {x) (12) 

which we iterate numerically. Here F[p(.),x] and po are 
defined in Eqs. (pX]); < u < 1 is a positive step pa- 
rameter. For convergence, the step parameter is usually 
set to a ~ 10 _1 . As pointed out in Rcf. S, the integral 
equation method outlined above allows to obtain density 
profiles corresponding only to thermodynamically stable 
or metastable states; for unstable states the iterations 
defined by Eq. ( p^ ) diverge for any a. 

The results in the form of s(e) vs. e and ds/de vs. e 
plots for various soft potential radii xq are presented in 
Figs jsf]l(i|. In Fig. || we show the entropy and inverse 
temperature plots for a relatively small soft potential ra- 
dius, xq = 10~ 2 . For this value of xq the system clearly 
exhibits all signs of gravitational phase transition. For 
comparison, in the same figure we present the /3(e) plot 
for a low-degeneracy fermionic system with p = 10 4 from 
Fig. 3. Despite the completely different nature of the 
short-range cutoffs for these systems, their uniform state 
entropies exhibit a strikingly similar behavior. This once 
again illustrates that properties of a uniform state are 
determined mostly by the long-range properties of the 
interaction. Naturally, the core-halo state structure and 
its properties depend on the nature of the cutoff, so the 
corresponding branches in the entropy vs. energy plot are 
visually different. Nevertheless, the difference is weak so 
the physical properties of phase transitions in long-range 
self- attracting systems are relatively insensitive to the 
precise form of the small-scale regularization. 

Similarly to the fermionic system, as the range of the 
cutoff is increased, the range of existence of metastable 



FIG. 8: Plots of entropy s(e) (top) and inverse temperature 
/3(e) = ds/de (bottom) vs. energy e for the soft potential 
radius xo = 10 -2 . Plot of /3(e) for the fermionic system with 
p = 10 4 is shown in dashed line. 
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FIG. 9: Plots of inverse temperature /3(e) = ds/de vs. energy 
e for different soft potential radii. 



states shrinks and finally disappears. At this point the 
gravitational phase transition crosses over to the first- 
order one. From the data presented in Fig. || we esti- 
mate that this crossover happens at xq — x gr ~ 0.021. 
For x gr < xq < x cr the system exhibits a normal first- 
order phase transition until the critical point (e cr ,x cr ) 
is reached. The plots presented in Fig. |l^ indicate that 
e cr w —0.7 and x cr ~ 0.22. When xq > x cri no phase 
transitions are present in the system. This allows us 
to conclude that similarly to the fermionic system, as 
the soft potential radius is increased, the self-attracting 
system with soft Coulomb interaction exhibits consec- 
utively gravitational, first-order, second-order, and no 
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FIG. 10: Plots ol inverse temperature /3(e) = ds/de vs. 
energy e for the soft potential radii (bottom to top) xo = 
0.12, 0.2, 0.22, and 0.25. 



phase transition at all. 



VALIDITY OF MEAN FIELD 
APPROXIMATION 



IV. 



To obtain the plots shown in Figs. [2|-|10| we used the 
saddle point or the mean-field approximation. It raises 
an important question of whether the distinct features 
of the gravitational phase transition described above are 
intrinsic or appear as artifacts of this approximation. Be- 
fore the approximation is applied, the microcanonical en- 
tropy S(E, N) is expressed through the logarithm of a 
sum of microscopic densities of states Wi of all macro- 
scopic states with the energy E and number of particles 
N (fcs = 1). The mean- field approximation is equivalent 
to replacing a sum of contributions from these macro- 
scopic states, usually represented by a functional integral 
over corresponding density (or phase space density) pro- 
files p, by a contribution from the single state or density 
profile po g 1,0: 



S = ln^2 Wi ~ In J Vp W[p] « \nW[p ] 



(13) 



This state po extremizes the density of states W and, 
consequently, the entropy S. The condition 



SW 
~Jp~ 



0. 



(14) 



defines po that gives a global maximum, a local max- 
imum, or a minimum (or saddle point) to the entropy 
which correspond to stable, metastable, or unstable 
states, respectively. Let us first consider the range of 



energies, (E < E 1 and E > E 2 in Figs. || and ||, or all 
energies in Figs. |?] and [l(], for which there exists only 
one global entropy maximum. This case is described 
in traditional textbooks of thermodynamics: the sin- 
gle equilibrium state corresponds to a very sharp max- 
imum in the density of states and the first-order cor- 
rections to the entropy per particle s scales as 1/N, 
i.e. s = s[po] + 0(1/N). When stable and metastable 
states coexist (E\ < E < E2), the mean- field ap- 
proximation is also an asymptotically exact approxima- 
tion for the stable states, since the relative contribu- 
tions from the metastable states p m to W[p] scales as 
exp{N(s[p m ] — s[po])} However, the sharp kink in the 
mean-field entropy plot at E* (as in Fig. l]) appears in 
the true, non-mean-field entropy plot only in the N — > 00 
limit; as for any finite number of particles the metastable 
states contribute significantly to W[p] in the vicinity of 
E* , where s[p m ] — s[po] ~~ * 0. Similarly, the mean-field 
approximation works well for the metastable states when 
they are sharp local maxima of W[p\. But this approx- 
imation breaks down when the contributions to the en- 
tropy from the metastable state p m becomes comparable 
to or less than the contributions from other states p' in a 
vicinity of p m , \\p m -p'\\ < II Pm II- This happens when p m 
ceases to be at least a local maximum of entropy, which is 
exactly what takes place at the metastability-instability 
transition points E± and E2. This breakdown of the mean 
field approximation near the E\ and E2 energies can also 
be viewed as fluctuation-induced uncertainty in the exact 
location of the metastability-instability transition. It is 
shown in Ref. |^0) that the relative uncertainty AE/Ei 
in the position of the collapse point E\ scales with the 
number of particles as TV -2 / 3 . Hence, given that N is 
large, the mean-field results are asymptotically exact for 
all energies except for those near the ends of metastable 
branches E\ and E 2 - 

Another distinct feature of gravitational phase tran- 
sitions is the anomalous stability of the metastable 
branches [E U E*\ and [E*,E 2 ] (Fig. Q). Consider, for 
example, a metastable uniform state and a stable core- 
halo state both having the same energy somewhere in the 
middle of the interval [E\,E*]. The entropy minimum, 
that separates the entropy maxima corresponding to the 
stable and metastable states, has the depth AS which is 
proportional to the number of particles, i.e. AS = NAs 
(for example, in Fig. || where As is equal to the differ- 
ence in coordinates between the metastable (solid line) 
and unstable (dashed line) states, As sa 0.1). Physically, 
this is so because the transition from a metastable uni- 
form state to a stable core-halo state requires a macro- 
scopic fluctuation equivalent to the rearrangement of the 
density distribution in the whole system. Hence, the 
probability of the metastable-stable transition is propor- 
tional to exp(— NAs) and becomes prohibitively small 
even for a moderate number of particles TV. Only near 
the ends of metastability branches E\ and E2 is the prob- 
ability of metastability-stability transition significant; it 
is of order O(N ) in the interval [E\,E\ + AE], where 
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uniform phase 



cr 



core-halo phase 



cutoff range 



FIG. 11: Sketch of the phase diagram of a self-attracting 
system. The crossover point between the gravitational and 
first-order phase transition and the critical point are marked 
"gr" and "cr", respectively. 



AE ~ E 1 TV- 2 / 3 ||. 

Given the arguments presented above we conclude that 
the mean-field approximation adequately represents the 
phenomenology of self-attracting systems and correctly 
describes the gravitational, first-order, and second-order 
phase transitions. 



CONCLUSION 



In Sect. |Tj and |TTl we considered two examples of 
self-attracting systems, the ensemble of self-gravitating 
fermions and the ensemble of classical particles interact- 
ing via attractive soft Coulomb potential. These sys- 
tems have a similar ~ 1/r interaction at large distances 
but very different forms of short-range cutoffs. While in 
the second example the short-range cutoff is evidently 
ro, in the first example the role of the short-range cut- 
off is indirectly played by the Pauli exclusion princi- 



ple, which depends on the particle density. Despite this 
different small-r behavior, both of the considered en- 
sembles exhibit the same sequence of phase transitions: 
gravitational, first-order, second-order, and none, as the 
range of their respective cutoffs is increased. The sketch 
of their phase diagram in cutoff-energy coordinates is 
represented in Fig. [n] and completes the one given in 
Ref. |l(| in cutoff-temperature coordinates. The ex- 
amples considered in this paper were chosen mainly be- 
cause of the physical importance of the 1 jt potential but 
similar phase diagrams exist for all nonintegrable 
< a < 3 attractive potentials. The main physical rea- 
son behind this phase diagram is that the short-range 
cutoff controls the maximum density a self-attracting sys- 
tem can achieve. As the range of the cutoff is increased, 
the collapsed central core becomes less dense and occu- 
pies more volume, and at some point the system has sim- 
ply nowhere to collapse. It happens when the central 
density of the non-collapsed uniform state becomes sim- 
ilar to the core density of the collapsed core-halo state. 
Likewise, the critical point is reached when the maxi- 
mum allowed density becomes so small that the system 
remains virtually uniform for any energy. The purely ge- 
ometrical nature of these arguments indicates that the 
phase diagram, obtained in the previous two sections, 
should be robust and insensitive to such simplifying as- 
sumptions as spherical symmetry. The validity of the 
main approximation used in this paper, the mean-field 
approach, was discussed in Sec. [V. It is revealed that 



the mean-field approximation correctly describes the be- 
havior of the self-attracting systems for all accessible en- 
ergies and cutoff radii, excluding the immediate vicinities 
of the collapse points. We leave the study of the collapse 
points in the finite-particle systems as well as the dynam- 
ics of the collapse for a future paper. 
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